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ABSTRACT 



The goal of this work is to study the evolution of high redshift (z 1.4) qui- 
escent galaxies over an effective area of ~ 1.7 deg 2 in the COSMOS field. Galaxies 
have been divided according to their star-formation activity and the evolution of the 
different populations, in particular of the quiescent galaxies, has been investigated in 
detail. We have studied an IRAC (mag^.etim <22.0) selected sample of ~ 18000 galax- 
ies at z ^ 1.4 in the COSMOS field with multi- wavelength coverage extending from 
the U band to the Spitzer 24 /im one. We have derived accurate photometric red- 
shifts (oa 2 /(i+ z ,) = 0.06) through a SED-fitting procedure. Other important physical 
parameters (masses, ages and star formation rates (SFR)) of the galaxies have been 
obtained using Maraston (2005) models. We have divided our sample into actively 
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star-forming, intermediate and quiescent galaxies depending on their specific star for- 
mation rate (SSFR = SFR/M). We have computed the galaxy stellar mass function 
(GSMF) of the total sample and the different populations at z = 1 .4 — 3.0. We have 
studied the properties of high redshift quiescent galaxies finding that they are old 
(f — 4 Gyr), massive ((M) ~ 10 10 ' 65 M Q ), weakly star forming stellar populations with 
low dust extinction (E(B-V) ^ 0.15) and small e-folding time scales (r ~ 0.1 — 0.3 
Gyr). We observe a significant evolution of the quiescent stellar mass function from 
2.5 < z < 3.0 to 1.4 < z < 1.6, increasing by ~ 1 dex in this redshift interval. 
We find that z ~ 1.5 is an epoch of transition of the GSMF: while the GSMF at 
z > 1.5 is dominated by the star-forming galaxies at all stellar masses, at z < 1.5 the 
contribution to the total GSMF of the quiescent galaxies is significant and becomes 
higher than that of the star-forming population for M ^ 10 10 75 M Q . The fraction of 
star-forming galaxies decreases from 60% to 20% from z ~ 2.5 — 3.0 to z ~ 1.4 — 1.6 
for M ~ 1O 11 '°M0, while the quiescent population increases from 10% to 50% at the 
same redshift and mass intervals. We compare the fraction of quiescent galaxies de- 
rived with that predicted by theoretical models and find that the Kitzbichlcr & White 
(2007) model, implemented on the Millennium Simulation, is the one that better re- 
produces the shape of the data. Finally, we calculate the stellar mass density of the 
star-forming and quiescent populations as a function of redshift and find that there 
is already a significant number of quiescent galaxies at z > 2.5 (p ~ 6.0 M Q Mpc~ 3 ), 
meaning that efficient star formation had to take place before that time. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: star formation 



1 INTRODUCTION 

Understanding the processes that regulate stellar mass 
growth in galaxies as well as tracing the history of galaxy 
star formation and mass assembly over the cosmic time 
jLillv et all Il996l. iMadau et all 1 19961 iPozzetti et"aH 1 19981 . 
iDickinson et all 120031 ') are among the most important and 
discussed topics in modern cosmology. During the last 
decade, extraordinary progress has been made in our com- 
prehension of the formation and evolution of galaxies. 
Thanks to deep multi-wavelength surveys we are now able 
to observe large samples of distant galaxies, thus reducing 
cosmic variance and allowing us to study objects at early 
cosmic times. 

In the present Universe, it is known that galaxies 
can be separated into two broad populations: passively 
evolving galaxies with red colors, also known as quiescent 
galaxies or early type galaxies (ETGs h ereafter), and star- 
forming galaxies with b lue colors (e.g.. [Baldrv et al]|2004l . 
iBrinchmann et all 20041). A similar division exte nds to at 
least z ~ 1 dBell et all 12004 IWillmer et all 120061) and pos- 
sibly to z ~ 2 jGiallongo et al.ll2005l . ICirasuolo et al.l [2007I . 
ICassata et alj|2008l ). 

ETGs are the most massive galaxies in the present day 
Universe and are characterized by simple and homogeneous 
properties (morphologies, colors, passively evolving stellar 
populations, scaling relations). As a consequence, they are 
a key population to investigate the stellar m ass assembly o f 
massive galaxies ove r the cosmic time (e.g.. [Renzinil [20061 ). 
iGallazzi et al.l |2006l ) studied the properties of local quies- 
cent galaxies and, in particular, the color-magnitude and 
the Mg2-cr„ relation. They found that for high-mass ellip- 
tical galaxies the dispersion in age is small, while at the 
low mass end there is a tail towards younger ages, reflect- 
ing a shift in stellar growth towards less massive galaxies 
in recent epochs. Besides, they also argue that at increasing 



stellar mass there is also an increase in both total metallicity 
and a/Fe ratio, suggesting that massive early type galaxies 
formed in a relatively short time-scale. The stellar mass and 
ages of the passive ETGs at z ~ 1—2 require precursors char- 
acterized by strong (> lOOMo/yr) and short-l ived (0.1 — 0.3 
Gyr) starbursts occurring at z > 2 - 3 (e.g., ICimattil [20091 
and references therein). 

Stellar mass assembly in galaxies subdivided by spectral 
and morphological types, as well as by star formation activ- 
ity has been inves t igated at z > by many authors (e.g., 
Bundv et al. Il2005l . iFranceschini et al. 120061 iPannella et al.l 



20061 . IPozzetti et alj|20ict ). ICimattil ^POd ) studied the B- 
band luminosity function of ETGs since z ~ 1, finding that 
the amount of the evolution for the ETG population de- 
pends critically on the range of masses considered, i.e., most 
of the massive ETGs are already in place at z ~1, while the 
density is still increasing with time for lower masses. This 
means that more massive galaxies have older stellar popu- 
lations and formed their stars earlier and more rapid l y tha n 
low mass galaxies, in agreement wi th IGallazzi et al.l ( 2006T ). 
Therefore, the downsizing pattern (|Cowie et al. 19961 ) could 
be applied to the assembly process of ellipticals at z < 1. The 
downsizing trend is also confirmed by a recent work based 
on a spectroscopic sa mple up to z ~ 1 in the COSMOS field 
|Pozzetti et alj|2010l ) that showed how the ETGs increase 
in number density with cosmic time faster for decreasing 
M. They also found that the number density of blue or spi- 
ral galaxies with M > 10 10 M© remains almost constant, 
while the most extreme population of star-forming galaxies 
at intermediate/high mass is rapidly decreasing in number 
density with cosmic time. The authors suggest a transfor- 
mation from blue active spiral galaxies of intermediate mass 
into blue quiescent and successively into red passive types 
with low specific star f ormation. Th i s is in agreement with 
the result obtained by iFaber et al.1 (|2007l ) who concluded 
that most present-day E/S0 near L* arose from blue galaxies 
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with ongoing star-formation that were quenched after z ~ 1 
and then migrated to the red sequence. The properties of 
nearby ETGs support a mixed scenario in which quenched 
galaxies enter the red sequence via wet, gas-rich mergers, 
followed by a limited number of dry, stellar mergers along 
the sequence. The evolution of the stellar mass function of 
the ETGs up to z ~ 1 has also been studied by other au- 
thors ( e.g-lBell et alj|2004 iBrown et ail l2007llTavlor et all 
l2009l , lNichol et al.ll2007l ~ finding a moderate increase in the 
number density of quiescent galaxies with cosmic time which 
is mass dependent. 

At higher reds hifts lArnouts et al.1 (|2007l ) and 
ICirasuolo et all (|2007f ) found a rapid rise in the space 
density of massive red sequence ga laxies from z ~ 2 t o z ~ 
1. These results are confirmed by lllbert et al.1 (|2010h who 
found that the stellar mass density of quiescent galaxies 
of all masses increases by 1.1 dex between z ~ 1.7 and 
z ~ 1, while it evolves only by 0.3 dex b etween 0.8 — 1 
and z ~ 0.1. This trend is also confirmed by I Font ana et all 
(|2009T l . who found that the fraction of quiescent galaxies 
increases from 15 — 20% at z > 2 to ~ 40% at z ~ 1.2, 
and recently bv lNicol et al.l (|201ll ). who found that the red 
sequence massive galaxies (M > 10 11 M©) increase in mass 
density by a factor ~ 4 from z ~ 2 to 1. 

Some groups have claimed the existence of a sub- 
stantial population o f massive galaxies at z > 4 (e.g., 
Mobasher et all 12004 lYan et al.1 [20061. IRodighiero et ail 
20071. iFontana et alj _ 120091 . llVIancini et_alj l2009fT 



Marchesini et all ( 20091 ) and iMarchesini et al ] l|2010l ) 



found that the number density of most massive galaxies 
(M > 3 x 10 11 Mq) seems to evolve very little from z ~ 4 
to z ~ 1.5, with a larger subsequent evolution down to 
z ~ 0.1. These results ar e broa dly consistent with the work 
from Ivan Dokkum et"afl (|2010h . where the authors derived 
an increase of massive galaxies by a factor of ~ 2 since 
z ~ 2 to the present. They concluded that this growth is 
likely dominated by mergers, as in situ star-formation can 
only account for ~ 20% of the mass build up from z — 2 to 

2 = 0. 

Several studies suggest that the critical redshift range 
where the stro ngest evolution and assembly took place is 1 
z s? 2 (eg.. lArnouts et al.l [20071 . lAbraham et al.1 120071 ). 
However, while the results up to z ~ 1 are quite solid 
thanks to the large samples of galaxies and to the spec- 
troscopic infomation available, the samples of ETGs studied 
up to now at z > 1.4 are still small. Besides, ETG are the 
more clustered galaxies in the universe, making their study 
very complicated given field to field variations which intro- 
duce large uncertainties due to the cosmic variance. Any- 
way, many works have identified these objects spectroscopi- 
cally, showing that passive/quiescent g alaxies with elliptical 
morphology may ex i st up to z ~ 2.5 (ICimatti et alj|2004 
iMcCarthv fc GDDSl |2004 ICimatti et alj|200gh . There are 



even some candidates of ETGs at z > 4 — 5, but their na- 
ture is still unknown as they are optically too faint for being 
spectroscopically observed. Therefore, the uncertainties in- 
volved in the study of the evolution of quiescent galaxies at 
high redshifts are still large due to the observational diffi- 
culty to identify large samples of ETGs at high z (> 1.5). In 
particular, long wavelength data are necessary to cover the 
optical-near infrared part of the spectrum of high redshift 
galaxies, a spectral range fundamental when calculating the 



galaxy stellar masses. In addition, the collection of a large 
sample of high redshift galaxies requires a combination of 
large areas and deep near-IR (NIR) observations. Due to 
these difficulties, the study of the galaxy stellar mass func- 
tion (GSMF hereafter) divided by galaxy types at z > 1.5 
has been explored only roughly and our knowledge of the 
evolution of galaxies at high z is still rather poor. 

Motivated by the limited information available at z ^ 
1.5, we present in this paper the evolution of the GSMF 
with cosmic time at z ^ 1.4 for the total galaxy popula- 
tion and for the sample divided on the basis of the galaxy 
star-forming activity in the COSMOS field. The COSMOS 
survey is the perfect wide area sample to look for rare 
and clustered objects such as the old and massive ETGs 
with much improved statistics. This survey (|Scoville et al.l 
120071 ) presents many advantages when compared with other 
surveys. First of all, the continuous coverage of ~ 2 deg 2 
substantially reduces the effect of cosmic variance. Besides, 
its extensive multi-A c overage allows to calculate accurate 
photometric redshifts jllbert et al.l 120091 ') and to study the 
sources based on their spectral energy distributions (SEDs), 
which is a more physical and complete approach than using 
only colors or morphological information. Secondly, the mul- 
tiwavelength data include deep Spitzer/I RAC (3.6 — 8.0/mi; 
Sand ers et al. 2007) and A" s -band data (|McCracken et all 
2010). crucial for estimating stellar masses at z ~2. 

Galaxies of the present work are selected in the IRAC1 
channel (3.6 (jm), allowing to sample the optical-NIR stel- 
lar bump at z > 1. In addition, data at red wavelength are 
crucial to detect galaxies that are very faint in the opti- 
cal (even undetected) but may contribute significantly to, 
or even dominate, the stellar mass density at high re dshift 
jYadl200d . fFranx et al.ll2003l . IRodighiero et al.ll2007fi . The 
extensive multi-band data are used to separate galaxies into 
different populations based on their SE Ds. The present wor k 
extends to hig her redshifts the st udy of lllbert et al.l (|2010T ) . 
The sample of lllbert et all (|2010h was selected at 3.6 /im as 
the present one, but their restriction to IRAC sources with 
an i + -band counterpart (i + < 25.5) limited their study to 
2 < 2 and to larger stellar masses. Our work was possible 
after a careful analysis of sources not detected in the optical. 

The paper is organized as follows: in Section 2 we intro- 
duce the COSMOS data and the cross-correlation method. 
In Sections 3 and 4 we describe the method to estimate the 
photometric redshifts and the galaxy stellar masses, while 
the criteria used to split the sample in different populations 
is described in Section 5. In Section 6 we study the main 
properties of the high redshift quiescent galaxies and in Sec- 
tion 7 we present the GSMF for the whole sample and for 
the different populations, discussing our results in Section 
8. Finally, in Section 9 we summarize our conclusions. 

Throughout this paper we use standard cosmology 
(fi m = 0.3, fi A = 0.7), with H = 70 km s" 1 Mpc" 1 . Magni- 
tudes are given in AB system. The stellar ma sses are given in 
units of solar masses (M ) for a Chabrier (|Chabrierl 120031 ) 
initial mass function (hereafter IMF). The stellar masses 
based on a Salpeter IMF were converted into a Chabrier 
IMF by adding —0.2 4 dex to the logarithm of stellar masses 
jArnouts et alj|2007t ) . 
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2 DATA 

We have built a multi-wavelength IRAC 3.6 (im selected 
catalogue in the COSMOS field by associating the IRAC 
sources with the optical and MIPS counterparts through 
the likelihood ratio method (see Section 2.2). For sources 
undetected at 5a in the i + -band selected catalogue we have 
performed a nearest match with the ^-selected catalogue. 
The IRAC catalogue is 95% complete at 5.75 fi3y (22.0 
mag). The final catalogue consists of 78649 sources with 
magzsum < 22.0: 95% (74742) of them with an optical 
counterpart, 4.5% (3554) with no optical counterpart but 
a T^s-band counterpart and 0.5% (353) with only IRAC de- 
tection. Besides, 14% (11352) of the IRAC sample are also 
detected in the MIPS band. 



2.1 Catalogues 

(i) IRAC 

The COSMOS IRAC catalogue (|Sanders et al.l I2007M in- 
cludes photometry in the 4 IRAC channels (3.6, 4.5, 5.8, 8.0 
/im ) for sources that have a measured flux in IRAC Channel 
1 above 1 fijy (23.9 mag) in the COSMOS field. We used 
the 1.9" aperture corrected flux and we considered sources 
with mog3.6 M m < 22.0. 

(ii) Optical 

The optical catalogue (|Capak et al.l 12007m is an i + -band 
(i < 26.5 mag) selected catalogue in the area of the Subaru 
survey (2 deg 2 ) and includes photometry in 9 bands (u* , Bj, 
g + , Vj,r + , i + , z + , J, K s ). We applied syste matic offsets 
derive d to optimise the photometric redshifts bv lllbert et al.l 
|2009l ) and auto-offset corrections (which account for the 
aperture correction) to the optical catalogue. 

(iii) 21 fxm 

The 24 /xm catalogue (|Le Floc'h et al.ll2009h is composed 
by sources detected above 0.15 mjy (18.5 mag). 

(iv) K, 

The 7f s -selected catalogue from iMcCracken et all (|2010l ) 
contains galaxies detected at K a and additional data from 
Bj,i + ,z + and J bands. The catalogue is complete down to 
K a ~ 23 mag . We applied a magnitude d ependent aperture 
correc tion to the fluxes, as explained in IMcCracken et all 
(|20ld ). 



2.2 The cross-correlation 

We have performed a cross-correlation of the IRAC sources 
with the optical and MIPS catalogues through the likelihood 

ratio method. 

The likelihood ratio technique (|Sutherland fc Saunders! 
Il992l ; ICiliegi et all 120031 ) has the advantage of taking into 
account not only the position (as for the nearest neighbour 
method, commonly used in literature), but also the mag- 
nitude distribution of the background sources, helping to 
reduce wrong identifications. The likelihood ratio LR is de- 
fined as the ratio between the probability that the source is 
the correct identification and the corresponding probability 
of a background, unrelated object: 



probability distribution function of the positional errors of 
the IRAC sources (assumed to be a two-dimensional Gaus- 
sian) and n(m) is the surface density of background objects 
with magnitude m. In order to derive an estimate for q(m) 
we have first counted all objects in the optical/MIPS cata- 
logue within a fixed radius around each source (total(m)). 
This has then been background subtracted and normalized 
to construct the distribution function of real identifications: 



q(m) = 



real(m)Q 
J2 t real(m), 



(2) 



LR = q(m)f(r)/n(m) 



(1) 



where q(m) is the expected probability distribution, as 
a function of magnitude, of the true counterparts, f(r) is the 



where real(m) = total(m) — n(m) (the sum is over all 
the magnitude bins of the distribution) and Q is the prob- 
ability that the optical/MIPS counterpart of the source is 
brighter than the magnitude limit of the optical/MIPS cat- 
alogue (Q = J" 1 '" 71 q(jn)drn). Finally, a best threshold value 
for LR (LR t h) is chosen in order to maximize both the com- 
pleteness (avoid missing real identifications) and reliability 
(keep spurious identifications as low as possible) of the as- 
sociations. 



2.2.1 Results 

In this Section we present the final multi-wavelength cata- 
logue created by cross-correlating the IRAC catalogue with 
the optical and the MIPS catalogues. For the sources with 
no 5cr i + -band counterpart we have performed a match with 
the J^a-band catalogue. 

(i) IRAC + Optical 

We considered sources in the optical catalogue with flux in 
the i + -band > 5a, and eliminated the masked areas and the 
stars from both catalogues. We performed a match through 
the likelihood ratio method considering a maximum radius 
of 3 arcsec (in order to maximize the overdensity around 
the IRAC position) and a dispersion of a = 1" for the dis- 
tribution of the positional errors. The distribution of the 
background, the total and the real counts is shown in Fig. 1. 

The majority (96%) of the associations have a separation 
$5 1 arcsec, with the median separation being ~ 0.34 arc- 
sec. For a non negligible fraction of the associations with 
a separation > 1 arcsec, the photometric redshift we ob- 
tained (see Section 3) was not in agreement with the spec- 
troscopic one, meaning a probable spurious IRAC-optical 
association. Therefore, we have considered as reliable only 
matches within a separation of 1 arcsec, which is a standard 
separation used in literature when cross-correlating optical 
and IRAC catalogues. To reduce contamination from Active 
Galactic Nuclei (AGN) we have also eliminated the IRAC 
sources with an optical counterpart asociated to an X-ray 
detection with XMM-Newton (1.5% of the IRAC sample). 
Finally, we obtained a reliable optical counterpart for 74742 
IRAC sources (95%). We will refer to these sources as the 
"optical sub-sample" hereafter. 

(ii) IRAC +24 fim 

We found 11352 IRAC sources (14%) with a reliable MIPS 
counterpart detected at 5a. The likelihood ratio method, as 
described above, has been used. 

(iii) IRAC+K S 

For the sources without an optical counterpart we have 
performed a match with the A^-selected catalogue. In this 
case we used the nearest neighbour technique with a maxi- 
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Figure 1. Magnitude distribution of the background (black line), 
total (green line) and real counts (red line) in the i+-band. 



mum separation of 1 arcsec. The similar wavelengths of the 
K s and IRAC1 bands allows to perform this match with 
a low probability of misidentification. 90% (3554 sources 
over 3907) of the IRAC sources with no optical coun- 
terpart, have a K 3 counterpart ( "K s sub-sample" here- 
after). As already ment ioned, the T^-band catalogue from 
iMcCracken et all (|2010t ) contains also information in addi- 
tional bands (Bj,i + ,z + and J). The median magnitude in 
the i + -band of these sources is i + ~ 26, i.e., these sources 
are optically very faint, consistent with them being not in 
the "optical sub-sample ". 

Finally, we have 353 sources for which we did not find 
neither an i + -band counterpart, nor a K 3 counterpart and 
for which we have only the IRAC bands information ( 
"IRAC sub-sample" hereafter) and, eventually, the 24 fim 
(50 sources). The results of the cross-correlation and the 
final numbers of our catalogue are summarized in Table 1. 



3 PHOTOMETRIC REDSHIFT 

The only efficient way to estimate the redshifts of the 
galaxies for the whole sample, given the large number of 






N (% of the whole sample) 


MIPS detection 


IRAC sample (all) 


78649 


11352 


Opt. sub-sample 


74742 (95%) 


10779 


K s sub-sample 


3554 (4.5%) 


523 , 


IRAC sub-sample 


353 (0.5%) 
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Figure 2. Representative templates of the SED library. Or- 
ange: simple ste llar population models of 0.05, 1.0 and 3.0 Gyr 
iMarastonll2005ri. Red: el liptical l lPolletta et alj2007T>. Green: spi- 
rals llPolletta et al .120071 ). Blue: blue galaxies |l3ruzual fc Charlotl 



200; 



i (F 
3). 



Table 1. Number of sources detected in the optical, K s and MIPS 
catalogues. 



sources considered, is to determine their photometric red- 
shifts through a detailed SED-fitting procedure. 

P hotometric redshifts were derived using the Le Phare 
code l|Arnouts et al. 1 12001| . lllbert et al.l 120061 ) , which com- 
putes photometric redshifts based on a simple x 2 template- 
fitting proce dure. We u s ed th e COSMOS SED library 
provided by lllbert et al. I (|2009l) . which contains a set of 
templates composed of ellipticals and spirals from the 
iPolletta et al.l (|2007l) library and also includes blue galax- 
ies from the Bruzual fc Charlotl (|2003h population synthesis 
models. We added 6 simple stellar population (SSP) models 
of differen t ages (0.05, 0.5, 1.0, 1.5, 2.0 and 3.0 Gyr) from 
iMarastorJ (|2005h to account for passively evolving galaxies 
at high redshifts. The final library is composed of 7 ellip- 
ticals, 12 spirals, 12 starbursts and 6 SSP templates. Some 
representative templates of the library are shown in Fig. 2. 
This set of templates maximizes the accuracy of the derived 
photometric redshifts (z p ) when compared with the available 
spectroscopic redshifts (z s ) (see Section 3.1.). 

W e allowed 3 different extinction laws |Calzetti et al.l 
|2000| , Calzetti modified and IPrevot et al] 1 19841 ) to be ap- 
plied to each template. The Calz etti modified exti nction law 
includes a UV bump at 2175 A filbert et al.ll2009l ). which is 
fundamental to have good z p at z > 1.5. We used 9 different 
values of E(B-V): 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4 and 
0.5. 

We set an upper limit (5<r) to each band when there is 
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K s sub-sample 
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22.0 
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1 -fl/lO 1 






IRAC2 


23.3 




IRAC3 


21.2 




IRAC4 


21.0 




MIPS 


18.5 





Table 2. 5<r limit adopted for each band 



no detection, as shown in Table 2. This 5a limit has been 
chosen following a conservative approach, i.e., as the bright- 
est magnitude value of the sources in the catalogue detected 
with S/N=5. We fitted our data from the U to the 24/im 
band. The \ 2 minimization procedu re is highly affect ed by 
the adopted photometric errors (see lllbert et al] 120091 for a 
detailed discussion of the issue). We have increased the flux 
error by 6, 5 and 8% for the optical, IRAC and MIPS bands 
respectively. 

3.1 Photometric Redshift accuracy 

We were able to test the qualitity of our derived photomet- 
ric redshifts for the IRAC sources with i + -band detection 
by comparison with a large sample of spectroscopic red- 
shifts. In particular, we made use of 8176 sources with de- 
rived spectroscopic redshift at a very high confidence level 
(99.8%, flags 3.1, 3. 5, 4.1, 4.5) from the z COSMOS faint 
and bright surveys (|Lillv et alil2007l , 120091 '). We estimated 
the redshift accuracy as o"a z /(i+z s ) (where Az = zp — zs 
is the difference between the photometric and the spec- 
troscopic redshifts) using the normalized median absolute 
deviation defined as 1.48(|z p — z a |/(l + z s )). We also de- 
fined the percentage of catastrophic errors, rj, as the objects 
with \z p — z s |/(l + z s ) > 0.15. We obtained an accuracy of 
CAz/(i+z B ) = 0.060 and a percentage of catastrophic fail- 
ures of 7? = 3 .3% a s can be seen in Fig. 3. As a comparison, 
lllbert et al] (|2009r ) using the bright (i + < 22.5 mag) sam- 
ple with 30 bands and narrow filters obtained a — 0.007 
and < 1% of catastrophic errors for their derived photo- 
metric redshifts, whil e the average value of a obtained by 
iPerez-Gonzalez et al. I l|2008h using an IRAC selected cata- 
logue with data from the UV to the MIR is 0.055. Given the 
fact that we use 14 broad bands, we find our result compa- 
rable to those in literature. 

We also tested the reliability of our z p for the faintest 
sources, with mag-ss^m > 21.0 and i + > 22.5, as shown in 
Figs. 4 and 5. The value of a for the faint i + -band sources 
(a = 0.057) is comparable to the value of a for the whole 
sample, while the number of catastrophic errors increases 
(rj = 12%). For the faint IRAC1 sample both a and rj are 
smaller than for the whole sample (a = 0.037, 77 = 1.9%). 
This may surprise the reader, but it must be noticed taht the 




Figure 3. Comparison between z p and z s for the 8176 sources 
with optical and IRAC detection of the zCOSMOS sample. The 
value of a and rj are shown in the plot. The different colors repre- 
sent sources fitted by different SED templates as explained in Fig. 
2. Dashed lines represent \z p — z a \/(l + z a ) = 0.15, i.e., sources 
out of that region are considered as catastrophic errors. 

values for the whole spectroscopic sample are highly affected 
by a significant dispersion on low redshift sources ( z s ~ 
0.3 — 0.5, see Fig. 3), which are present in a smaller fraction 
in the faint IRAC sample. 

3.2 Photometric Redshift accuracy for the IRAC 
sources with no I band counterparts 

For the sources with no i + -band counterpart it is not possi- 
ble to directly compare the photometric redshifts with spec- 
troscopic redshifts, since no z s are available for these sources. 
To test the quality of the derived photometric redshifts we 
have made use of the la dispersion of the z p distribution, 
given as an output of the Le Phare code. The mean error 
on the z p for the whole sample is a z = 0.17, which corre- 
sponds to <J&. z /n+ Za ) ~ 0.057 at z ~ 2, meaning that the 
la dispersion of the z p distribution is consistent with the 
°"Az/(i+z s ) derived from the z p and z s comparison. The red- 
shift uncertainties derived from Le Phare have been used in 
the Montecarlo simulations (see Section 7). 

In Figs. 6,7 and 8 we show some examples of SED fit- 
ting for sources fitted by different templates (elliptical (Ell), 
simple stellar populations (SSP), spirals (S) and blue galax- 
ies or starbursts (SB)) for sources from the optical, K s and 
IRAC sub-samples. 

3.3 Redshift distribution 

Fig. 9 shows the redshift distribution for each sub-sample, 
highlighting the contribution from different classes of best- 
fitting templates. Templates have been divided in 3 classes: 
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Figure 4. Comparison between z p and z a for 183 sources with 
i+-band > 22.5. Colors and dashed lines as in Fig. 2. 




12 3 4 



Z phot 

Figure 5. Comparison between z p and z s for 2381 sources with 
mags.Q^m > 21.0. Colors and dashed lines as in Fig. 2. 

elliptical/old SSP (older than 1 Gyr), blue/young SSP 
(younger than 1 Gyr) and spirals. It can be observed that 
these redshift distributions are not smooth, but show some 
peaks at given redshifts, e.g., at z ~ 1.0, 2.0 for the K s 
sub-sample or a very prominent peak at z ~ 4.0 for the 
IRAC sub-sample. Note that the peak at z ~ 4 in the 
IRAC sub-sample is not used in this paper since we limit 
our analysis to z < 3.0 and also the small number of sources 



ident 238462 ident 238398 




0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 

X[/im] X[*im] 



Figure 6. Example of SED-fitting for 4 sources from the optical 
sub-sample. The observed SED of each source (black filled cir- 
cles) is shown with the corresponding best-fit solution (red solid 
line), as well as the derived z v and best template model. Arrows 
represent upper limits (5a) when no detection is found in that 
band. 



ident 168697 ident 161984 




0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 

fyml Mem] 



Figure 7. Example of SED-fitting for 4 sources from the K 3 
sub-sample, i.e., s ources not present in the i + -band catalogue 
llllbert et al,ll201oT ) but with a K 3 counterpart fMcCrack cn et all 
l20ld) . 

involved in this peak (< 0.2% of the sample ). This kind 
of features in the redshift distribution are presumably not 
real, but due to the determination of the redshifts through a 
SED fitting procedure. Spiky redshift distributions have also 
been repor ted in literature, e .g., by IPerez-Gonzalez et alj 
(|2008h and iFranx et ail (120081 ) among others. The number 
of sources in each subsample with a reliable z v is shown in 
Table 3, as well as the number of sources for each sub-sample 
fitted with each class of template. The number of high red- 
shift (z > 1.4 and z > 2.5) sources for each sub-sample is 
also shown. 

It is clear how differently the three sub-samples behave, 
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Figure 8. Example of SED-fitting for 4 sources from the IRAC 
sub-sample, i.e., sources only detected in the IRAC bands. 



OPT sub-sample 



OPT sub-sample 



OPT sub-somple 




Figure 9. Rcdshift distribution for the different sub-samples. In 
the 3 upper panels we show the rcdshift distribution of the optical 
sub-sample (black line), as well as the redshift distribution of the 
sources fitted by an elliptical/old SSP template (red line), a spiral 
template (green line) or a blue/young SSP template (blue line). 
The same is plotted for the K s and the IRAC sub-samples in the 
middle and bottom panels, respectively. 



with 88% and 91% of the K s and the IRAC sub-samples re- 
spectively found at high redshift [z ^ 1.4), in comparison 
to 21% of the optical sub-sample. For very high redshift 
sources (z > 2.5) the differences are even more significant, 
since only 2.3 % of the optical sub-sample is found at those 
redshifts, while this percentage is 27% and 47% for the K s 
and IRAC sub-samples, respectively. Also the relative frac- 
tions of best-fitting templates are very different. On the one 





Opt. sub-sample 


K s sub-sample 


IRAC sub-sample 


All 


74742 


3554 


328 


Ell/old SSP 


10% 


12% 


22% 


o 


25/0 


4o/o 


71 07 


Blue/young SSP 


62% 


43% 


7% 


z > 1.4 


21% 


88% 


91% 


z > 2.5 


2.3% 


27% 


47% 



Table 3. Percentages of best fit templates and high redshift 
sources (z ^ 1.4 and z 2.5) for the different sub-samples. 



hand, 10% of the optical sub-sample is fitted by an ellip- 
tical/old SSP model, while for the IRAC sub-sample the 
fraction is double (22%). On the other hand, only 7% of 
the IRAC sub-sample is fitted by a blue/young SSP model, 
while these templates represent 62% of the solutions for the 
optical sub-sample. These percentages are indicative of the 
importance of an IRAC selected sample to study the prop- 
erties of high redshift galaxies whose SED is reproduced by 
quiescent populations. 



4 GALAXY STELLAR MASSES 

To study the evolution of the galaxy stellar mass function 
of sources at high redshift we have selected the sources with 
z ^ 1.4. When available, we have made use of z s (654 out 
of 19042 sources, ~ 3.4% of the sample). 
The galaxy stellar masses have been derived by means of the 
Lephare code by fitti ng our data (up to 5.8^m) with a set of 
SED templates from IMarastonl (|2005l ) with star formation 
histories exponentially declining with time as SFR tx e~ t//r . 
We used 9 different values of r (0.1, 0.3, 1.0, 2.0, 3.0, 5.0, 
10.0, 15.0 and 30.0 Gyr) with 221 steps in age. The metal- 
licity is solar an d the IMF Chabrier. Dust extinction was 
applied using the lCalzetti et all (|2000l ) extinction law, with 
a maximum E(B-V) value of 0.5. We imposed to the derived 
age of the galaxies to be less than the age of the Universe at 
that redshift and greater than 10 8 years (the latter require- 
ment avoids having galaxies with extremely high specific 
star formation rates SSFR = SFR/M (SSFR hereafter). 
The IMarastonl (|2005l ) models include a better treatment of 
the thermally pulsing asymptotic giant branch (TP-AGB) 
phase, which has a high impact on modelling the templates 
at ages in the range 0.3 < t < 2 Gyr, where the fuel con- 
sumption in this phase is maximu m, specially for the near- 
IR part. Although some authors dKriek et al.l |2010l) have 
recently claimed that the iMaraston ( 20051 ) models do not 
properly reproduce the observed data in a sample of post- 
starburst galaxies, i.e., at the time when TP-AGB stars are 
thought to be mo st dominant, w e have instead found good 
fits for our SEDs. iPozzetti et all (|2010l ) measured an aver- 
age systematic shift of 0.14 dex betwee n the s tellar masses 
co mputed with the iBruzual fc Charlotl (|2003l ) models and 
the IMarastonl (|2005l ). We will take into account this differ- 
enc e when comparing our res ults with those obtained with 
the lBruzual fc Charlotl (|2003l ) models. 

Uncertainties in the stellar mass derivation are due to 
a number of different assumptions in the SED fitting proce- 
dure, like, for example, the use of different IMF, extinction 
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laws, metallicities, star formation histories (SFH) or SED 
template libraries. Tests on simulated catalogues consider- 
ing the effect on stellar mass estimates of different choices of 
reddening law, SFHs, metallicities and SED libraries show 
a typical dispersion o f the order of a(logM) ~ 0.20 (see 
iBolzonella etal .11 20091 ), showing that the stellar mass is a 
rather stable parameter in SED-fitting. 

The mean error on the stellar mass, as estimated from 
the la dispersion of the mass distribution (given as an out- 
put of the Le Phare code), is 0.16 dex, consistent with the 
uncertainties mentioned above. 



5 GALAXY CLASSIFICATION 

We are interested in studying the characteristic evolution 
of the quiescent and star-forming galaxies. To discriminate 
between these two types, we have made use of a classification 
based on the SSFR derived through the SED-fitting process. 
We considered 3 different populations: 

(i) Actively star-forming galaxies: those for which 
logiSSFRlGyr- 1 ]) > -0.5 

(ii) Intermediate galaxies: those for which —2 < 
logiSSFRlGyr- 1 ]) < -0.5. 

(iii) Quiescent galaxies: those for which 
log(SSFR[Gyr- 1 ]) < -2. 

We found 84 sources (~ 0.5% of the high redshift sample, 
i.e., sources with 1.4 < z < 3.0) with log^SSFRlGyr' 1 ]) < 
—2 and detected in the 24/im band. As the emission at 24 /im 
usually arises from dust heated by star-formation activity, 
we include these sources in the intermediate population. 

The limits on the SSFR values have been defined from 
the distribution of the SSFR values (Fig. 10 for the opti- 
cal sub-sample) and we have verified that they are consis- 
tent with the best-fitting template classification. Most of 
the sources with very low SSFR are fitted by an elliptical 
or old SSP template, while those with the highest activity 
are generally fitted by a young SSP or a blue template, as 
should be expected on the basis of their colors. Considering 
the whole sample, we have carefully studied the cases where 
star-forming sources are fitted by an elliptical or old SSP 
template, as well as quiescent sources fitted by a blue or a 
young SSP template. We found only 26 sources (< 0.1% of 
the high redshift sample) fitted by an elliptical or old SSP 
template with high activity. On the other hand, more than 
94% of the 639 sources (3.5% of the high redshift sample) 
fitted by a blue template with low SSFR are not detected in 
the U band, which is a critical band to determine the star 
formation activity of a galaxy. 

In Table 4 we report the number of galaxies classified 
as star- forming, intermediate and quiescent, as well as the 
MIPS detected sources, in each redshift bin for sources with 
1.4 < z sC 3.0. 

Other tests, such as the BzK classification (|Daddi et al.l 

I2004T ). which separates star forming and quiescent galaxies 
with respect to their (B — z), (z — K) colors for sources with 
1.4 ^ z ^ 2.5, shows that galaxies for which low activity is 
derived roughly occupy the passive region in the BzK dia- 
gram ((z - K) - (B - z) < -0.2 and (z - K) > 2.5), while 
for the star-forming galaxies, their colors are consistent with 
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all 


quiescent 


intermediate 


star-forming 


MIPS 


1.4-1.6 


5142 
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1325 
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4346 
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2.0-2.5 


4693 


287 


914 


3492 


1178 


2.5-3.0 


1794 


77 


415 


1302 


272 



Table 4. Number of galaxies with 1.4 ^ z s: 3.0 in each redshift 
bin as classified by the SSFR. 
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Figure 10. SSFR distribution for the optical sub-sample of 
sources with z ^ 1.4. The colors represent the different classes 
derived from the SED fitting templates. Dashed vertical lines rep- 
resent the limits considered in SSFR to divide sources as actively 
star-forming, intermediate or quiescent. 

the BzK star-forming selection (Fig. 11). Only 3 sources de- 
fined as active star forming fall in the passive region of the 
BzK diagram. 

The lack of detection in some bands worsens the derived 
SSFR, thus complicating a precise classification. We have 
derived mean SSFR errors making use of the la dispersion 
of the SSFR distribution, obtaining a mean SSFR error of 
0.38 dex for the whole sample. 



6 PROPERTIES OF HIGH-Z QUIESCENT 
GALAXIES 

The evolution over the cosmic time of quiescent galaxies rep- 
resents a key tool to understand the galaxy mass assembly. 
In this Section we summarize the properties of the high red- 
shift quiescent galaxies (also known as ETGs in literature) 
by analyzing the physical parameters obtained through the 
SED-fitting procedure. 

We find 1455 quiescent galaxies at z ^ 1.4. They are old 
and massive galaxies ((M) ~ 10 10 ' 65 Mq): more than 56% of 
them have M ^ 1O 1O ' 6 M and ~ 12% have M ^ 1O 11O M . 
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Figure 11. Bxif diagram llDaddi et alj|2004h for optically de- 
tected sources with 1.4 ^ z ^ 2.5. The colors represent the differ- 
ent classes derived from the SSFR criteria (red: quiescent, green: 
intermediate, blue: star- forming). The dashed lines show the dif- 
ferent regions of the BzK diagram. 







ident 20135 








ident 23252 




100.00 


.Moss=11.09 






100.00 


.Mass=10.35 






10.00 

i 

| 1.00 


E(B-V) = 0.00 
S5FR = -2.58i' 
.ogc = 3.00C* 






10.00 

1 

2 i.oo 


E(B-V) = 0.00 
SSFR = -2.20 






0.10 


T = 0.3GyrJ 


OPT 




0.10 




OPT 




0.01 








0.01 








a 


1 1.0 


10.0 


100.0 


0.1 1 


10.0 


100.0 






ident 16364 








ident 33312 




100.00 


.Mass= 10.95 






100.0 


.Mass = 1 1 .23 






10.00 

I 

0.10 
0.01 


E(B-V) = 0.10 
SSFR = -2.69 

r = 0.1Gyr jT 


K 




10.0 

1 

0.1 


E(B-V) = 0.00 
5SFR = -4.56 


K 


2 = 2.628. 





1 1.0 


10.0 


100.0 


0.1 1 


10.0 
X|>n] 


100.0 




Mass= 10.63 


ident 66687 








ident 240468 




10.00 


E(B-V) = 0.15 
SSFR = -3.23 




z = 2.567. 


10.00 


E(B-V) = 0.15 
SSFR = -6.98 






I ,.00 


.nqe=1.?8Gyr 








.age=2.20Gyr 






0.10 


T = 0.1Gyr 1 > 


IRAC 




0.10 


T = 0.1Gyr 1 


/ IRAC 




0,01 








0.01 








a 


1 1.0 


10.0 


100.0 


0.1 1 


10.0 


100.0 



Figure 12. SED-fitting of some high redshift quiescent galaxies. 
Also shown the main physical parameters derived through the 
SED-fitting process (mass, extinction, SSFR, age, t, z). The OPT, 
K s and IRAC labels indicate the sub-sample at which each source 
belongs. 




Figure 13. Stellar mass as a function of redshift (small grey 
dots). Red crosses are the Mn m of the 20% faintest galaxies at 
each redshift, while M m i„(z) are plotted as yellow circles. 



They have weakly star-forming stellar populations, with 
ages ranging from 1 to 4.25 Gyr. Some of them are so old to 
reach the limit allowed by the age of the Universe at their 
redshift, meaning very high formation redshifts. They have 
e-folding time scales r ~ 0.1 — 0.3 Gyr, i.e., they should have 
formed in a very intense and brief starburst. They are char- 
acterized by very low dust extinction (E(B-V) ~ — 0.15), 
implying that their red optical colors are due to evolved stel- 
lar populations and not to dust extinction. These c haracter- 
istics are in agreement with previous results (e.g., ICimattil 
l2009l . lKriek et al.ll200Sh . 

In Fig. 12 we show the SED-fitting and main physical 
parameters of 6 different high redshift quiescent galaxies. We 
show examples from the optical, K s and IRAC sub-samples. 
We underline the importance of the upper limits in the op- 
tical and K s bands for the IRAC sub-sample to constrain 
the SED model. 



7 THE GALAXY STELLAR MASS 
FUNCTIONS 

We study the evolution of the quiescent galaxies by deriving 
their mass function in four redshift bins. The mass function 
of the whole sample and of the star forming and intermediate 
galaxies are also estimated to compare the evolution of the 
different p opulations. Th e classical non-parametric 1/Vmax 
formalism (|SchmidtJll968l ) has been used and data have been 
fitted to a Schechter function (|Schechterlll976t l. The galaxy 
stellar mass function GSMF is computed only in the non- 
masked regions with a total covered area of 1.73 deg 2 . 

Estimating the stellar mass limit for a magnitude lim- 
ited sample is not straightforward due to the high range of 
possible M/L ratios for different galaxy populations and col- 
ors. To account for this effect we define, at each redshift, a 
minimum mass, M m i„, above which the derived GSMF is 
essentially complete because all types of galaxies are poten- 
tially observa ble above this mass. 

Following [Pozzetti et al] (|2009T ). we calculated for each 
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galaxy its limiting stellar mass, i.e., the mass it would 
have at its redshift if its observed magnitude were equal 
to the limiting magnitude of the sample (mag-i.uim, = 22.0); 
log(M Um ) = log(M) + 0A(mag 3 . 6 - mag ZMim ). In order to 
derive a representative limit for our sample we used the Mu m 
of the 20% faintest galaxies at each redshift. We then defined 
the M m in (z) as the upper envelope of the Mu m distribution 
below which lie 95% of the Mu m values at each redshift. 
This M m i n corresponds to a 95% completeness limit at each 
redshift. The distribution of mass with redshift, as well as 
the Mum and the M m i„(z) calculated as explained above 
are shown in Fig. 13 . 

Fig. 14 presents the total and the quiescent GSMF 
(QSMF hereafter) of our IRAC selected sample in 4 different 
redshift bins: 1.4-1.6, 1.6-2.0, 2.0-2.5 and 2.5-3.0. At z > 3.0 
the uncertainties become dominant so we decided to limit 
our analysis to z ^ 3.0. The circles represent the GSMF for 
all the sources, while the lines show the best Schechter func- 
tion fitting our data. The Poissonian errors have been c om- 
puted assuming Poissonian statistics |Zucca et etl.ll 199?! ^ : 



However, the Poisson errors in our GSMF are an under- 
estimate of the real uncertainties, since they do not take into 
account how the z p uncertainties propagate on the GSMF. 
This uncertainty is especially important for the optically un- 
detected sources, therefore we have to be very careful in the 
z-bins where these sources dominate (higher redshift). To as- 
sess for this uncertainty we have performed Montecarlo sim- 
ulations. We have created 20 mock catalogues by randomly 
picking a redshift within the redshift probability distribution 
function (PDFz) of each object (given by Le Phare code), 
we have recalculated the mass and other important physical 
parameters and we have computed the GSMF for each of the 
20 mock catalogues. We have then calculated the la disper- 
sion of the GSMF values at each mass and redshift bin. The 
cosmic variance is another important source of un certainty. 
Following the formalism of ISomerville et all (|2004h we have 
derived a cv ~ 0.18,0.14,0.14 and 0.16 respectively in the 
four redshift bins. Finally, we added in quadrature these er- 
rors to the Poissonian errors (yellow area in Fig. 14). The 
red area represents the upper and lower limits of the QSMF 
calculated in the same way as for the total GSMF. 

7.1 Comparison of the total GSMF with literature 

In Fig. 14 we show for comparison the GSMF 
derivations from th e litera ture (blue triangles from 
Perez-Gonzalez et alj |2008l , green asterisks from 
Fontana et al.1 12006I . red diamonds from iMarchesini et al.l 
20091 . orange squares from lllbert et alj|2010h . All GSMF's 
were converted to Chabrier IMF and shifted by —0.14 dex 
in mass (to account for the d iffer ence in the derived mass 
when using iMarastonl 120051 or iBruzual fc Charlotl 120031 
models, see IPozzetti et al.l 201dl ). We are in reasonably 
good agreement with many of the other results. However, a 
significant difference is observed in all redshift bins at high 
masses with respect to the results of Perez-Gonzalez et al. 
(2008), as they found a nu mber density sign i ficant ly higher 
than ours. With respect to IMarchesini et al.l (|2009l ). we find 
less massive galaxies at z < 2 and more intermediate mass 



2 

'— ' -5 

2 




1.4<z<1.6 

1.6<z<2.0 
2.0<z<2.5 
2.5<z<3.0 



10.0 10.5 11.0 11.5 12.' 

logM [M ] 

Figure 15. Evolution of the QSMF with redshift for different 
redshift bins. Blue: 1.4 — 1.6; green: 1.6 — 2.0 ; orange: 2.0 — 2.5; 
red: 2.5 — 3.0. Shaded regions represent the upper and lower limits 
of the GSMF for each redshift bin (see Section 7), while the thick 
lines are the Schechter fitting to the data. 



(logM ~ 10 . 5) gala xies at 2 < z < 2.5. The agreement with 
lllbert et al.l l)2010h . whose sample is also IRAC selected 
in the COSMOS field, is very good, with small differences 
probably arising from different red s hift binning, as it 
may happen also with iFontana" et al.l (|2006T ). However, we 
have a slight overdensity of intermediate mass galaxies 
(M ~ 1O 1O - 7 M ) at 2.0 < 2 < 2.5 when compared with all 
other results from literature, which may be too big to be 
explained only by differences in redshift binning. 

In comparison with other authors, we are using data 
from the COSMOS survey, which is the largest studied area 
(~ 1.7 deg 2 ) with full multi-wavelength coverage. 



7.2 The evolution of the quiescent galaxies 

In Fig. 15 the QSMF obtained in the four studied redshift 
bins is reported. The shaded areas represent the upper and 
lower limits, derived as explained in Section 7, while the 
thick lines are the Schechter fitting to our data. 

We observe a significant evolution of the QSMF from 

2.5 < 2 < 3.0 to 1.4 < 2 < 1.6, increasing ~ 1 dex for 
galaxies with logM ~ 11.0 in the redshift interval studied. 
The evolution is ~ 0.3 dex in each redshift bin, meaning 
that the number of quiescent galaxies continously increases 
with cosmic time. 

In Fig. 16 we compare the GSMF derived for quies- 
cent galaxies, with those of the other populations. The red, 
green and blue areas represent the quiescent, intermediate 
and actively star-forming populations, respecti vely. We also 
show for comparison the results obtained by lllbert et al.l 
with purple, dark green and light blue circles rep- 
resenting their quiescent, intermediate and high activity 
galaxies respectively. The agreement is good for the quies- 
cent population in both redshift bins in common (1.4 — 1.6, 

1.6 — 2.0), while the differences for the intermediate and 
star-forming population can be explained by the different 
classification methods used. While for the quiescent popu- 
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Figure 14. Stellar mass function in 4 redshift intervals for the quiescent (red area) and the total (yellow area) sample. Our estimates at 
each redshift for the whole sample are plotted as filled black circles. Data are fitted with a Schechter function (solid lines). The colored 
areas represent upper and lower limit of the total and quiescent mass functions as derived from the combination of the Montecarlo 
simulations, cosmic variance and Poissonian errors. Black dotted line: local MF iCole et al.ll20oJ) . Dash ed vertical lines represent th e 
completeness in mass for each redshift bin . GSMFs from other authors are also plotted: blue triangles from lPerez-Gonzalez et_al.l 
green astrerisks from lFontana et al, ] ll2006h , red diamonds from iMarchesini et al.l ||2009T ). orange squares form lllbert et al.l d201of i. 



lation we used a similar cut in SSFR, for the cut between 
the intermediate and high activity populations they used 
log(SSFR[Gyr~l]) ~ —0.3, while we are using a value of 
-0.5. 

The evolution of the star-forming galaxies is quite com- 
plex (see Fig. 16). They remain almost constant from 2.5 < 
2 < 3.0 to 2.0 < z < 2.5, being the dominant population 
at all masses in those redshift bins. However, their num- 
ber density decreases when moving to the lowest redshift 
bins and this decrease depends on the galaxy mass: while 
for high masses we observe a significant decline (~ 0.8 dex 
for logM ~ 11.0), for low masses (logM ~ 10.5) the GSMF 
is almost constant. A decrease of the star-forming popula- 
tion from 1.6 < z < 2.0 to 1.4 < z < 1.6 is in agreement 
with the results of lllbert et al.l (|2010t ) a nd, as reported in 
studies at lower red shift (see for example lllbert et al.ll201ol . 
iPozzetti etHI 120091 ). the number of star-forming galaxies 
continues to decrease to almost local redshift (z ~ 0.1). The 



intermediate population constantly increases its number at 
all masses, first at high masses (~ 0.2 dex from 2.5 < z < 3.0 
to 1.6 < z < 2.0 for logM ~ 11.0), then at low masses (~ 0.2 
dex from 1.6 < z < 2.0 to 1.4 < z < 1.6 for logM ~ 10.5). 
This means that they become more abundant than the star- 
forming population at logM > 10.6 at 1.4 < z < 1.6. Finally, 
as mentioned before, for the quiescent population we find an 
increase of ~ ldex in the number density from 2.5 < z < 3.0 
to 1.4 < z < 1.6 for logM ~ 11.0. At the lowest redshift bin 
( 1.4 < z < 1.6) the quiescent population becomes more 
important than the active one for masses logM > 10.75. 

Considering all the populations, we find that z ~ 1.5 is 
a clear epoch of transition of the GSMF: while the GSMF 
at z > 1.5 is dominated by the actively star-forming galax- 
ies at all stellar masses, at z < 1.5 the contribution to the 
total GSMF of the quiescent galaxies is significant and both 
the intermediate and the quiescent population become more 
important than the star-forming one for M ^ 1O 1O ' 75 M0. 
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Figure 16. Evolution of the GSMF for different activity types galaxies. The colored areas represent the upper and lower limits (see 
Section 7) while the thick lines are the Schechter fitting to the data. Red: quiescent. Green: intermediate. Blue: actively star-forming. 
The circles r eprese nt the GSMF for quiescent (purple), intermediate (dark green) and high activity galaxies (light blue), as obtained by 
lllbert et all ll2010h . Also shown arc the limits in mass (vertical dashed lines). 



In Table 5 we report the best-fit Schechter parameters 
for the quiescent sample. In the last redshift bin (2.5 < z < 
3.0 ) the Schechter fit failed to converge when leaving all the 
parameters free. We fixed the faint end slope a to its value 
at 2.0 < z < 2.5 to be able to obtain the other Schechter 
parameters. While, given the large uncertainties, we cannot 
derive firm conclusions on the faint end slope a nor on M* , 
we detect, as already mentioned, a significant increase of (f> 
from 2.5 < z < 3.0 to 1.4 < z < 1.6. 

We caution the reader that the uncertainties of this de- 
rived GSMF at such high redshift may change our results 
significantly due to the difficulties in estimating the z p as 
well as the activity classification. 



8 DISCUSSION 

To better understand the behaviour of the different popula- 
tions we have plotted the fraction of quiescent/star-forming 
galaxies over the total versus mass for different redshift bins 



(Fig. 17). The fraction of the actively star-forming galax- 
ies decreases from ~ 60% to ~ 20% from 2.5 < z < 3.0 
to 1.4 < 2 < 1.6 for M > 1O 11O M , while that of quiescent 
galaxies increases from ~ 10% to ~ 50% in the same redshift 
and mass intervals. The evolution is also mass-dependent, 
with the most massive galaxies becoming quiescent first, 
while at M < 1O 1O,5 M0 the star- forming population still 
dominates the sample at all redshifts. 

This trend suggests that actively star-forming galax- 
ies, which dominate the sample at high redshifts (mainly at 
z > 2.0), may quench their star formation and become less 
abundant at later time. This quenching of the star forma- 
tion seems to be mass dependent, with the most massive 
galaxies quenching their star-formation first, while the less 
massive galaxies continue to form stars until later times. If 
the evolution of the star-forming galaxies should be that of 
becoming intermediate galaxies to finally end their lives as 
quiescent galaxies, the increase of this latter population at 
low redshifts could be explained, with the quiescent popu- 
lation being the dominant population at low redshift, as all 
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Type z-bin a log(M*)[M e ] </>[lCT 4 Mpc~ 3 dcx" 1 ] 

Quiescent 1.4-1.6 -0.08± 0.98 10.57 ±0.25 1.9± 0.5 

1.6-2.0 -0.41± 0.92 10.61 ±0.19 1.0± 0.4 

2.0-2.5 -0.34± 4.83 10.50 ±0.84 0.9± 0.8 

2.5-3.0 -0.34 (Fixed) 10.52 ±0.18 0.5± 0.6 



Table 5. Schechter parameters of the QSMF at 1.4 sg z ^ 3.0. 




10.2 10.4 10.6 10.8 11.0 11.2 11.4 



logM [M ] 



Figure 17. Upper panel: Fraction of quiescent galaxies over total 
for different mass values. Colors represent the different redshift 
bins as in Fig. 15. Lower panel: Fraction of the star-forming pop- 
ulation over total versus mass. Colored vertical dashed lines are 
the mass limits for each redshift bin. 



galaxies may evolve to that type once their SF has stopped. 
Therefore, we could be unveling the formation epoch of the 
very active galaxies and their following evolution from the 
blue cloud to the red sequence. 

Fig. 18 shows the evolution of the fraction of quies- 
cent galaxies with redshift for sources with M > 10 10 ' 85 M Q 
(red circles). The errors are a combination of the Monte- 
carlo simulations and Poissonian errors. For comparison, we 
also plotted the values obtained by Fontana et al. (2009, 
orange triangles) and from Damen et al. (2010, light blue 
dia monds), as well as the predictions of theoretical mod- 
els dMenci et alj|2006l. M06 here after; IKitzbichler fc White! 
120071 . K07; lNagamine et al.ll2006l , N06; [Fontanot et aL^OCm 
F07). We observe a good agreement between the fraction 
of quiescent g a laxies derived in this work and those from 
iFontana et afl (|2009T ). The values of iDamen et al] (I2010TI 
are al so consistent with our results, although Damen et al.l 
(|2010T l sources are selected with M > 1O 11O M and are de- 
fined as quiescent when SSFR < l/(3t.tr). The fraction of 
massive quiescent galaxies increases from ~ 8% to ~ 35% 
from 2.5 < z < 3.0 to 1.4 < z < 1.6, with the main evolution 
occuring between 1.6 < z < 2.0 and 1.4 < z < 1.6 (where 
the fraction increases from ~ 18% to ~ 35%). This con- 
firms the expected cosmological increase in number density 
of massive quiescent galaxies with cosmic time. All the the- 
oretical models agree in predicting a gradual increase with 
time in the fraction of galaxies with low SFR, although large 



Figure 18. Fraction of quiescent galaxies (log(SSFR[Gyr 1 ]) < 
-2) with M ~S> 10 10 - 85 M Q as a function of redshift. Red cir- 
cles are the results obtained in this work (errors are the com- 
bination of the Poisson errors and th e Montecarlo s i mulat ions), 
orange triangles are the v alues from IFontana et al ] l|2009h and 
light blue diamonds from IDamen et al, I tepid) . Lines refer to 
the predictions of theoretic al models as decribed in the legend: 
Menci et all ||2006|) CM 06'), IKitzbichler fc White! (l2007h (K07), 
Nagamine et al] l|2006r j (N061 and lFontanot et al] ^20071 ) (F07). 



discrepancies between data and theoretical models are ob- 
served. It is beyond the scope of the present paper to give 
a detailed explanation of the differences between models, 
but we will discuss some generalities. Some models under- 
predict the fraction of quiescent galaxies at all redshifts, 
such as the purely semi-analytical models (M06, F07), and 
in particular predict virtually no such objects at z > 2, in 
contrast to what is observed. This lack of massive galax- 
ies at high redshift for the mode l of F07 has already been 
noticed by iFontanot et al] (|2007T ) , who admitted that the 
downsizing trend of the galaxies is not fully reproduced yet 
in these models, and suggested that some kind of feedback 
mechanism could help to reproduce the data. A slightly 
higher fraction of quiescent massive galaxies is predicted by 
M06, as they include AGN feedback mechanism, due to the 
growth of supermassive black holes and the AGN triggered 
by interactions in the host galaxies. This AGN feedback en- 
hances the fraction of galaxies populating the red branch 
of the color distribution and is particularly effective at high 
z. However, it is not effective enough to reproduce the ob- 
ser ved data at z < 2.0. T he purely hydro-dynamical model 
by iNagamine et al] (|2006h . represented with a shaded area 
for three different timescales r of the star-formation rate 
(ranging from 2 x 10 7 yrs to 2 x 10 8 yrs), appears to over- 
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Figure 19. Evolution of the stellar mass density of quiescent 
galaxies. The result from this work are the red symbols. The val- 
ues represented by the squares have been calculated as the inte- 
gral of the Schcchtcr function fitting the data, while the diamonds 
are a mere sum of the data for all masses. Orange symbols are 
results from literature as explained in the article . We also show 
for comparison the evolution of the stellar mass density of the 
star-forming plus the intermediate galaxies (blue shaded area). 
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Figure 20. Evolution of the stellar mass density of star-forming 
plus quiescent galaxies. The result from this work are the dark 
blue symbols. The values represented by the squares have been 
calculated as the integral of the Schechter function fitting the 
data, while the diamonds are a mere sum of the data for all 
masses. Light blue symbols are results from literature as explained 
in the article. We also show for comparison the evolution of the 
stellar mass density of the quiescent galaxies (red shaded area) . 



predict the fraction of high mass quiescent galaxies. The 
disagreement in the stellar mass density at z > 1 between 
these mod els and the observatio ns has already been com- 
mented bv lNagamine et al.l|200(jft . The semi-analytical ren- 
dit ion of the Millennium N-b ody dark matter simulation 
by licit zbichler fc White! l|2007l ) is the only one which agrees 
with our observed data, being able to reproduce both the 
normalization and the shape of the increase in the frac- 
tion of quiescent massive galaxies with decreasing redshift. 
This model introduces radio mode feedback from the central 
galaxies of groups and clusters, which seems to be fundamen- 
tal for the predictions to be consistent with the observations. 

Finally, the non negligible fraction of quiescent galax- 
ies at z > 1.8 implies that these galaxies assembled most of 
their stellar mass either during an active star-burst phase 
or through important merging processes at higher redshifts. 
These star formation episodes must be quenched either by 
efficient feedback mechanism and/or by the stochastic na- 
ture of the hierarchical merging process. 

In Fig. 19 and 20 we have calculated the evolution of 
the stellar mass density with cosmic time for the quiescent 
sample (red symbols) and the intermediate plus the star- 
forming samples (dark blue) . The values have been obtained 
in two different ways. The data plotted as squares represent 
the integral of the Schechter function fitting our data, in- 
cluding the extrapolations to the fainter masses. The dia- 
monds are the sum of the data at all masses. As expected, 
the latter values are lower (since data are incomplete at low 
masses) but the two estimates are consistent within ~ 0.18 
dex. We do not show the value of the stellar mass den- 
sity for the intermediate and star-forming galaxies in the 
redshift bins 2.0 — 2.5 and 2.5 — 3.0 derived through the 
integration of the Schechter parameters, neither the same 
value at 2.5 < z < 3.0 for the quiescent galaxies, due to 



the high uncertainties in the derived parameters given the 
incompletness at low masses. In the plots we also show a 
colored area representing the upper and lower limits of the 
quiescent (red) and star-forming (blue) populations, as well 
as various results from literature. Orange symbols repre- 
sent the stellar mass density of the quiescent galaxies as 
derived by other authors, while light blu e symbols stand 
for t he star-formin g ones (squa r es fro m iKochanek et alj 
200 li diamonds f rom lDriver et ail (l2006f), plus symbo ls from 
Bell et al.) (120031). triangles from I Arnouts et al .1 20071. aster- 



isks fr om Franceschini et all 2009. cr osses from [Borch et all 
(|2006T) . dots from lllbert et alj|2010h . We are in very good 
agreement with Ilbert et al. (2010) in the common redshift 
bins. We observe that the increase in stellar mass density 
in the redshift range of interest (t ~ 2 — 4.5 Gyr) is very 
rapid for both populations (~ 0.5 dex for the star-forming 
plus intermediate populations and ~ 0.9 dex for the quies- 
cent one), then slowing down at lower redshifts. However, 
this slow down happens earlier in time for the star-forming 
population ( z ~ 1.2, t ~ 5.0 Gyr) than for the quiescent 
galaxies, which continue to rapidly assembly mass until later 
times [z ~ 1.0, t ~ 6.0 Gyr ), reaching the star-forming 
galaxies mass density at z < 1.0. 

The existence of a non negligible population (logp\M.Q 



Mpc~ 



6.0 at z ~ 2.7) of quiescent galaxies with high 



masses which have already undergone major star formation 
even at the highest redshifts (z > 2.5) is fundamental for our 
understanding of the galaxy formation processes and crucial 
for testing theoretical scenarios. 



9 SUMMARY AND CONCLUSIONS 

We derived the GSMF and stellar mass density in a 1.73 deg 2 
area in the COSMOS field from z = 1.4 to z = 3.0 for the 
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whole sample and for the sample divided into populations, 
paying special attention to the quiescent one. The GSMF 
estimate is based on ~ 18000 galaxies with mag^.e^m <22.0 
and z p ^ 1.4. We summarize our results below: 

(i) We find that an IRAC selected sample is fundamental 
for studying quiescent galaxies at high redshift. About 88% 
and 91% of the K s and IRAC sub-samples are found at z > 

1.4 respectively, with this fraction dropping to 21% for an i + - 
band selected sample. Considering the galaxy classification, 
62% of the optically selected sub-sample is fitted by a blue 
template, while the same fraction is much lower (~ 7%) for 
the IRAC sub-sample. 

(ii) We study the main properties of the high redshift 
quiescent sample (z ^ 1-4), finding that they are old and 
massive galaxies, with (M) ~ 10 10 ' 65 Mq and ages ranging 
from ~ 1 to ~ 4 Gyr. They have small e-folding time scales 
r ~ 0.1 — 0.3 Gyr and very low dust extinction (E(B-V) 
~ — 0.15), meaning passively evolving populations. 

(iii) We observe a significant evolution of the QSMF from 

2.5 < z < 3.0 to 1.4 < z < 1.6, amounting to ~ 1 dex 
for galaxies with logM ~ 11.0. The evolution is ~ 0.3 dex 
in each redshift bin, meaning that the number of quiescent 
galaxies continously increases with cosmic time. 

(iv) We find that z ~ 1.5 is a clear epoch of transition 
of the GSMF: while the GSMF at z > 1.5 is dominated by 
the star-forming galaxies at all stellar masses, at z < 1.5 the 
contribution to the total GSMF of the quiescent galaxies is 
significant and the quiescent galaxies become more impor- 
tant than the star- forming population for M ^ 1O 1O ' 75 M0. 

(v) The fraction of star-forming galaxies decreases from 
60% to 20% from 2.5 < z < 3.0 to 1.4 < z < 1.6 for 
M ~ 1O 11O M0, while the quiescent population increases 
from 10% to 50% at the same z and mass intervals. 

(vi) All of the theoretical models agree in predicting a 
gradual increase with cosmic time in the fraction of galaxies 
with a low SFR. However, only the K07 Millennium-based 
model is able to properly reproduce the shape of the data. 

(vii) We find a significant number of quiescent galaxies 
already in place at z > 2.5 ( p ~ 6.0 MqMpc -3 ), meaning 
that these galaxies assembled most of their stellar mass in 
previous epochs in an active star-burst phase or through 
important merging processes at higher redshifts. 
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